Spectroscopy of three strongly coupled flux qubits 
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We have carried out spectroscopic measurements of a system of three strongly coupled four- 
junction flux qubits. The samples studied cover a wide range of parameters with the coupling energy 
between neighboring qubits varying between 0.75 GHz and 6.05 GHz. The observed complicated 
spectra agree well with eight-level theory. The experiments are relevant for the realization of a 
tunable coupling between qubits. 
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The potential realization of a full scale quantum com- 
puter requires the ability of coupling multiple qubits to- 
gether preferably so that the coupling can be turned 
on and off at will. In the context of Josephson junc- 
tion qubits there is a number of promising theoretical 
suggestions^'^ as well as already several experiments 
with coupled qubita^'^'SiSifiiiLii. In order for a quan- 
tum computer to be truly scalable it must be possible to 
couple many qubits together without degrading the co- 
herence time severely. Tunable coupling has not been so 
far demonstrated in an experiment where the coherence 
time would be equally good as in the case of so-called op- 
timally biased qubitsiSiiii^ii^. A problem common with 
many coupling methods is that in order to realize a two- 
qubit gate the biases need to be switched away from the 
region in the parameter space where the decoherence is 
minimal. In other words, while the coupling between two 
(or more) qubits is strengthened, the coupling between 
the qubits and the environment is, in many cases, also 
strengthened. 

This paper describes spectroscopic experiments on 
three strongly coupled flux qubits, which is relevant for 
instance for the scheme suggested in Ref. In that 
scheme the parametric coupling of two detuned optimally 
biased flux qubits is realized through the microwave mod- 
ulation of their tunable mutual inductance realized using 
a third qubit. Thus, the system we study is a set of 
three antiferromagnetically coupled flux qubits. When 
the flux threading the qubits is near half flux quantum 
(modulo flux quantum $o) the system of three qubits is 
reasonably well described by the Hamiltonian 
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where ej = 2Ipj{^j — $o/2) is the energy bias of qubit 
j controllable through the flux threading the qubit 
loop. Near the half-flux-quantum point each qubit ex- 
periences a double well potential and the tunneling en- 
ergy through the potential barrier separating the wells 
is Aj. The wells correspond to currents of magnitude 
Ipj circulating in opposite directions along the loop and 
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FIG. 1: SEM image of the sample and the principle of the 
measurement. 



the above Hamiltonian is actually written in this circu- 
lating current basis. The antiferromagnetic interaction 
between the qubits k and / is characterized by the cou- 
pling strength Jki = Mkilpklpi where Mki is the mutual 
inductance. 

A sample realizing such a system is shown in Fig. ^ 
The interesting part of the sample consists of three four- 
Josephson-junction flux qubits, similar to those studied 
individually in Refs.H^andllsl coupled together by shar- 
ing one edge of the superconducting loop. The qubits 1 
and 2 are thus expected to be weakly coupled due to small 
geometric inductance whereas each of them is expected 
to be strongly coupled to the qubit 3 through kinetic in- 
ductance. Each qubit has one smaller Josephson junction 
whose area is about a « 0.5 times smaller than that of 
the larger junctions having areas of about 200 nm by 400 
nm. The qubits arc coupled to a four-junction readout 
SQUID consisting of the three larger loops. The SQUID 
can be seen as a dc SQUID whose junctions are replaced 
by dc SQUIDs. This design was chosen in order to re- 
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duce the influence of the inevitable noise in the bias cur- 
rent. The couphng of the qubits to the SQUID is also 
through kinetic inductance. Since the different states of 
the qubits correspond to different magnetic field configu- 
rations in the SQUID loops we expect a finite population 
of the excited states to result in a small but detectable 
change of the switching current of the SQUID. 

The samples used in this study were fabricated out 
of Al using standard shadow evaporation through a Ge 
enforced mask patterned with e-beam lithography. The 
designs of samples A, B and C were slightly different. 
In sample A we used an on-chip capacitor to shunt the 
SQUIDiiiSii^. In samples B and C this capacitor was 
removed in order to clean up the measurable spectrum 
from resonances such as LC resonances and the plasma 
resonance. An additional benefit of removing the large 
superconducting shunt capacitor was improved flux sta- 
bility. Both designs had no on-chip bias resistors but 
instead we used an on-chip LC filter consisting of a long 
superconducting lineii and a superconducting capacitor. 
The two plates of the Al parallel plate capacitor were de- 
fined in a separate e-beam lithography step and the insu- 
lator was formed by heavily oxidizing the bottom layer. 
The estimated cutoff of the filter was 100-200 MHz. In 
sample A the inductor was about 15 mm long and 600 
nm wide resulting in a transmission line resonance around 
10 GHz. The shorter line of samples B and C (3.2 mm) 
gives a resonance around 20 GHz. The oxidation of the 
junction was done using a mixture of O2 (10%) and Ar 
(90%) between the depositions of the 20 nm and 30 nm 
Al layers. For the microwave line and bonding pads we 
used evaporated gold film defined by optical lithography. 
Part of the microwave line also had an Al layer. The 
measurements were carried out in a dilution refrigerator 
with a base temperature around 20 mK. The fridge had 
a triple /i-metal magnetic shield. 

The principle of the measuremen1ii& is shown in Fig. ^ 
First a microwave pulse of typically 5 /iS duration is ap- 
plied to the microwave line on chip. Then immediately 
after this a current pulse of about 5-10 ns is applied to 
the SQUID. The height is chosen so that the SQUID 
switches about every second time. The current pulse 
is followed by a trailing plateau whose height is about 
70% of the switching current. The application of cur- 
rent through the SQUID changes the magnetic fiuxes 
experienced by the qubits and this shift is assumed to 
be adiabatic. The purpose of the plateau is to main- 
tain the SQUID in the voltage state if and only if the 
SQUID switched. By repeating each measurement typi- 
cally 10^ times and counting the relatively slow voltage 
pulses (few ^s) we can deduce the switching probability 
Psw under particular circumstances. To supply the fiux 
bias to the qubits we used an external coil capable of in- 
ducing about 20 Gs. To carry out spectroscopy on the 
qubits we first located them in the flux space. The ba- 
sic measurement was then carried out by sweeping the 
microwave frequency and magnetic field and repeating 
the above mentioned measurement scheme. At each fiux 
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FIG. 2: Spectroscopy of sample A. Top panel is an inten- 
sity plot of Psw while the bottom panel includes the nu- 
merically calculated spectrum. Black dots denote excita- 
tion energies from the ground state |0) for the states I*]/) 
with K0|X;j=iO-i|*)P > 0.002. In the calculation we used 
Ai/h = 2.2 GHz, Aa/ft = 2.2 GHz, As/h = 2.5 GHz, 
Jn/h = J23/h = 2.05 GHz, J12 = GHz, A-l-i/c&o = 
-0.0128, A$2/$o = -0.0102, Jpi = 194 nA, /p2 = 194 uA, 
^p3 = 179 nA. The lines shared by the qubits have the di- 
mensions of 1.35 ^mx30 nmxl60 nm. Qubits 1 and 2 have 
a = 0.5 while qubit 3 has a — 0.475. Oxidation was done 
with 25 mTorr for 5 min. Switching measurement yields 

point the height of the current pulse was adjusted to get 
roughly the same Pgw in the absence of the microwave. 
The qubits (as well as other resonances) cause deviations 
from this probability enabling the characterization of the 
spectrum of the three-qubit system. 

FigureElis an example of such a spectroscopic measure- 
ment on sample A. The horizontal lines are resonances 
due to the presence of on-chip capacitors. The strong res- 
onance around 10 GHz is most likely a half- wavelength 
resonance in the long bias line. However, the qubits are 
clearly visible and form a rich spectrum. The top panel 
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FIG. 3: Spectroscopy of sample B. The parameters used 
in the calculation are Ai//i — 1.3 GHz, A2//1 = 1.0 GHz, 
Az/h = 3.6 GHz, Ji3/h = J2s/h = 6.05 GHz, J12 = 
GHz, A$i/$o - -0.0096, A$2/$o = -0.010, /pi = 156 nA, 
Ip2 = 156 nA, /p3 = 124 nA. We only show levels for which 
!(0|Ej=i'^il*)P > 0.0001. The lines shared by the qubits 
have the dimensions of 1.35 /imx20 nmxlOO nm. Qubits 1 
and 2 have a = 0.5 while qubit 3 has a = 0.45. Oxidation was 
done with 25 mTorr for 10 min. We get /|S,ai° ~ 20.2 fiA. 



shows the measured spectrum only while in the bottom 
panel a theoretical calculation of the spectrum is shown 
on top of the measured spectrum. The calculation is 
simply done numerically using Eq. ^ by finding the ex- 
cited state eigenvalues and subtracting the ground state 
eigenvalue from them. It is worthwhile stressing that the 
theoretical part is not based on a fitting procedure such 
that there may be considerable amount of error in all 
the parameters, but the agreement does seem at least 
qualitatively very good. The persistent currents denoted 
by Ipk in the figure captions are the actual values used 
in the computation and do not account for the shield- 
ing effect. The actual persistent currents are expected 
to be significantly larger and also the mutual inductance 
is expected to be smaller than what it seems consider- 
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FIG. 4: Spectroscopy of sample C. In the calculation we 
used Ai/h = 5.0 GHz, Ai/h = 7.5 GHz, As/h = 8.2 GHz, 
Ji3/h = 0.79 GHz, Jvi/h = 0.755 GHz, J12 = GHz, 
A-l-i/c&o = -0.0041, A<I'2/$o = -0.0087, Ipi = 81 nA, /p2 = 
78 nA, /p3 = 64 nA. Levels with |{0| E^=i ff^l*>P > 0.005 are 
shown. The different slopes are due to the fact that the qubits 
1 and 2 were designed with a = 0.5 and a = 0.475 while the 
qubit 3 has a — 0.45. The lines shared by the qubits have the 
dimensions of 1.35 /imx20 nmxlOO nm. Oxidation was done 
with 35 mTorr for 10 min. We get I^^J,° « 5AfiA. 



ing the currents Ipk- Based on the measured normal 
state resistivity of 2K = 11 /iilcm we estimate usingi^ 
^kin = ($oei?„)/(7r2ABcs) that M13 = M23 = 31 pH. 
In order for this to be consistent with the coupling en- 
ergy the shielding should be about 11%. Due to the cou- 
pling we observe more than three levels and everywhere 
the spectrum lines cannot be simply associated with a 
particular qubit. Roughly speaking, the line having a 
minimum of about 2.5 GHz corresponds to the qubit 3 
and the other two minima around 2.2 GHz correspond 
to the qubits 1 and 2. This can be justified by not- 
ing that the qubit 3 is the only one coupled strongly to 
two other qubits. The slightly higher tunneling energy 
is because the qubit 3 was fabricated with a = 0.475 as 
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the designed ratio of the large and small junction area 
while qubits 1 and 2 had a — 0.5. Since we are using 
only one control magnet we can write the flux biases as 
<i>j = $ + A$j. The small offsets A<f>j are, in the ab- 
sence of trapped vortices, due to slight differences in the 
areas. Since we only work in rather small range around 
<I>o /2 the effect of area difference on the slope is not sig- 
nificant. We choose a convention where A<i>3 = 0. Note 
however that due to the interaction between qubits, not 
even the qubit 3 has minimum exactly at <&o/2. The lack 
of visibility of some levels at certain field can be under- 
stood as small transition matrix element. In fact, in the 
theoretical calculation we have plotted only levels I^E") for 

which K0|Ej=iO-^|4')P > 0.002. 

Figures 13 and 0] illustrate measured spectra for sam- 
ples B and C. As can be seen the calculations again 
agree extremely well with the measurement. Together 
with sample A the samples cover a very wide range of 
parameters. It is especially noteworthy that the cou- 
pling energy can be made as large as 6 GHz without 
using coupling junctions^. The shielding effect as well 
as the coupling inductance are expected to be very large 
in these samples since the line shared by the qubits and 
SQUID is made thinner and narrower than in sample A. 
A prediction based on normal state resistance yields the 
estimate M13 = M23 = 74 pH. In order to explain the ob- 
served coupling energy the shielding in sample B should 
be about 40% and in sample C about 14%. These num- 
bers are reasonable although hard to verify. The scaling 
of these percentages is correct as a function of max^- 

Even though the predicted kinetic inductances are 
large one may suspect that the large coupling could 
be partially due to shared Josephson inductance Lj. 
Namely, a fraction of the persistent current of each qubit 
could flow through the large junctions of the SQUID. 
However, this effect may only reduce the coupling from 
the prediction based on shared kinetic inductance since 
the effect of the large shared Lj is cancelled by the 
fact that the fraction of the persistent currents flow- 



ing through the big junctions decreases with increasing 
Lj. Explicitily, if the currents of two qubits k and I can 
flow through two inductances 2Li (kinetic) and 2L2 (ki- 
netic plus Josephson), and if half of these inductances 
are shared with the other qubit, then Jki = LiL2/{Li + 
L2)lpklpi- This clearly is maximal for L2 —>■ 00. An ad- 
ditional experimental point related to this is that in all 
samples we reproduced the measured spectrum well by 
setting the direct coupling of the qubits 1 and 2 to zero 
even though their currents have a chance to flow through 
a "shared" large Josephson junction. 

In conclusion we have carried out spectroscopic mea- 
surements of three strongly coupled flux qubits. We 
demonstrated that one can achieve a wide range of cou- 
pling strengths using kinetic inductance for the coupling. 
The relatively high coupling energy is attributable to the 
rather large normal state resistivity of our aluminum film. 
For instance in Ref. the resistivity seems to be about 
7 times smaller while e.g. in Ref. it is about half of our 
value. We find that the effective three-qubit Hamiltonian 
describes the measured spectrum well. Based on our ex- 
perience it is possible to reproducibly fabricate qubits 
whose flux biases differ by less than 1% even when us- 
ing only a single bias coil. The different effective area of 
qubit 3 in all samples is due to the fact that its position is 
not symmetric with respect to other qubits. However, it 
is possible to tune the area very accurately by adjusting 
the layout. The similar areas enable optimal point bias- 
ing of many qubits without strongly coupled independent 
bias lines. The experiments reported in this paper are a 
step towards the tunable coupling of flux qubits. The 
parameters of sample C in particular are very close to 
those required for the parametric coupling using a third 
qubit^. However, for good enough coherence the cou- 
pling to the readout will have to be modified to be more 
symmetric than in the present design. 

We would like to thank M. Grajcar for useful discus- 
sions. 
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